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Plasmons, the collective oscillations of interacting electrons, possess emergent properties that 
dramatically alter the optical response of metals. We predict the existence of a new class of plasmons 
- chiral Berry plasmons (CBPs) - for a wide range of two-dimensional metallic systems including 
gapped Dirac materials. As we show, in these materials the interplay between Berry curvature and 
electron-electron interactions yields ehiral plasmonie modes at zero magnetic field. The CBP modes 
are confined to system boundaries, even in the absence of topological edge states, with chirality 
manifested in split energy dispersions for oppositely directed plasmon waves. We unveil a rich CBP 
phenomenology and propose setups for realizing them, including in anomalous Hall metals and 
optically-pumped 2D Dirac materials. Realization of CBPs will offer a new paradigm for magnetic 
field-free, sub-wavelength optical non-reciprocity, in the mid IR-THz range, with tunable splittings 
as large as tens of THz, as well as sensitive all-optical diagnostics of topological bands. 


In electronic systems, chirality expresses the system’s 
ability to discriminate between forward and backwards 
propagation of electronic signals along certain directions. 
This technologically useful and hotly sought-after prop¬ 
erty can be achieved through the application of exter¬ 
nal magnetic fields. However, the need for strong ap¬ 
plied fields “on-chip” brings many challenges for appli¬ 
cations. Recently, materials exhibiting chirality in the 
absence of a magnetic field have started to gain promi¬ 
nence. These include metals exhibiting anomalous- [T] 
and quantum anomalous- nn Hall effects, as well as 
non-magnetic materials pushed out-of-equilibrium, where 
for example a zero-field charge Hall effect was recently 
demonstrated [7]. In each case, zero-field chirality arises 
from Bloch band Berry curvature, a fundamental prop¬ 
erty of Bloch eigenstates that dramatically affects single¬ 
particle electronic motion and material responses miHiEi. 

Here we show that Berry curvature can work in con¬ 
cert with interactions, leading to new types of collec¬ 
tive modes in two-dimensional (2D) “topological” metals, 
with non-vanishing Berry flux (i.e., net Berry curvature), 
T. In particular, T gives rise to chiral plasmonie excita¬ 
tions — propagating charge density waves with split dis¬ 
persion for oppositely directed modes — in the absence 
of a magnetic field (Fig. [^. We refer to these collec¬ 
tive modes as “chiral Berry plasmons” (CBPs). Notably, 
these chiral modes are localized to the edge of the 2D 
metal, even in the absence of topological edge states, and 
exhibit a rich phenomenology (see below). 

We expect CBPs to be manifested in a wide variety of 
magnetic as well as non-magnetic materials. The former 
are materials that exhibit anomalous Hall effects, wherein 
time reversal symmetry (TRS) breaking is encoded in the 
Bloch band Berry curvature. The latter include a range 
of readily available gapped Dirac materials, e.g., transi¬ 
tion metal dichalcogenides and graphene/hBN, wherein 
TRS breaking is achieved by inducing a non-equilibrium 
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FIG. 1: Chiral Berry plasmons (CBPs) at zero mag¬ 
netic field, a) CBPs are manifested along the edges of a 2D 
metal with non-vanishing Berry flux, T (inset). Counterprop- 
agating modes exhibit a split dispersion, with the 

splitting increasing with wave vector q along the edge. Above 
a threshold frequency cuth (shaded region) the fast mode 
merges with bulk plasmon modes, see Eq. ^ and 

dashed line, leaving a single unidirectional mode prop¬ 

agating along the edge, b) Zoom-in of (a) showing separa¬ 
tion of and below uth. c) and d) Electric po¬ 

tential for bulk and edge plasmon modes. Eor bulk modes 
the anomalous velocity field for electrons (indicated by ar¬ 
rowheads) runs parallel to the wave fronts and does not affect 
the speed of collective propagation. Near edges the build-up 
of surface charges leads to an anomalous velocity flow that 
assists the collec tive propagation for the mode, see text 

below Eq. (S-2). Eor the anomalous vel ocity flow op¬ 

poses the collective motion, see text below Eq. ( |S-2| . Param¬ 
eter values used: 7^ = 1.0 and ao = 3 for (a) and (b), see 
Eq. §, and \q/qQ\ — 0.3, T — 0.3, and ao = 3 for (c) and (d). 


valley polarization [7]. In both cases CBPs are charac¬ 
terized by clear optical signatures such as split peaks in 
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optical absorption. 

From a technological perspective, CBPs in non¬ 
magnetic materials are particularly appealing since they 
provide an entirely new platform for achieving a range of 
magneto-optical effect analogues that are magnetic field- 
free, and “on demand.” A prime example is optical non¬ 
reciprocity uni, which is central for optical device com¬ 
ponents, e.g., optical isolators and circulators. Above 
a threshold frequency (cjth^ see shaded area in Fig. [^), 
the single unidirectionally propagating mode allows 
for chiral transport of light via coupling to CBPs. Such 
CBP mediated waveguides provide a novel paradigm for 
deep sub wavelength [IIHI3], linear, and magnetic field 
free strong non-reciprocity, crucial for miniaturizing op¬ 
tical components. In particular, we predict that CBPs 
can enable non-reciprocity over a large technologically 
important bandwidth (THz to mid-IR). 

The intrinsic chirality of CBPs starkly contrasts with 
that achieved via cyclotron motion of charged particles 
in a magnetic field, with important quantitative as well 
as qualitative consequences. In the latter case, chiral¬ 
ity arises via the Lorentz force and gives rise to conven¬ 
tional magnetoplasmons inHn]. There, the cyclotron 
frequency, fiujc = TieBjm [18], determines the constant 
splitting between magnetoplasmon modes of opposite 
chirality, which can be of order a few meV for accessible 
field strengths. In contrast, chirality in CBPs arises from 
the combined action of plasmonic self-generated electric 
fields and the anomalous velocity of Bloch electrons, the 
phase space dual to the Lorentz force. This combination 
makes the CBP mode splitting directly sensitive to plas- 
mon wavelength, the Berry flux, and interaction strength, 
in contrast to magnetoplasmon splittings which only de¬ 
pend on magnetic field [To] [43] . As a result, for short 
wavelengths and unscreened interactions, large splittings 
hAuj of several tens to a hundred meV can be achieved 
[see Eq. (§ below]. 

Self-Induced Anomalous Velocity - The origin of 
CBPs in two dimensions can be understood from the Eu¬ 
ler equations for electron density, n(r, t) [20] : 

dtn{r, t) -f V • v(r, t) = 0, (la) 

dtp(r, t) - n(r, t) eV0(r, t) = 0, (lb) 

where v(r, t) and p(r, t) are the velocity and momentum 
density fields, (/)(r) is the scalar electric potential, and 
—e < 0 is the electron charge. We note that in principle 
the force equation &) also includes contributions arising 
from the stress density of the electronic fluid. However, 
at long-wavelengths these contributions yield only sub¬ 
leading corrections to the plasmon dispersion [21]. As 
our aim here is to clearly and most simply demonstrate 
the existence and main features of CBPs, in this work we 
neglect small corrections due to the stress density. Eor an 
alternative formulation in terms of currents and conduc¬ 
tivities, see supplementary online information (SOI) [22] . 


In order to fully specify the dynamics, Eq. 0 must be 
supplemented by a set of constitutive relations which re¬ 
late velocity, momentum, density, and the electric poten¬ 
tial. Plasmons emerge from Eq. 0 as self-sustained col¬ 
lective oscillations of n(r,t), with the potential 0(r, t) = 
J d‘^r' W{r^r')6n{r'^t) generated by the plasmon’s den¬ 
sity fluctuations 6n = n(r, t) —uq. Here no is the average 
carrier density, and IF(r, r') is the Coulomb interaction. 

As we argue, in the presence of Berry curvature, the 
constitutive relations take on an anomalous character. 
This can be seen by starting with the quasiparticle semi- 
classical equations of motion [8], Vp = ^ + ^p x r2(p), 
p = eV0(r), where Vp = f is the quasiparticle velocity 
and 5p and = f)(p)z are the Bloch band disper¬ 

sion and Berry curvature, respectively [44]. The cor¬ 
responding velocity density fields are found from these 
relations and the phase space density /^(r, p,t) by sum¬ 
ming over all momentum p and bands {i}, (5(r, t) = 
f /{2Trhr. This gives: 
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( 2 ) 


where T = / d‘^pCli{p)f^{p)/{27rh)^ is the (dimen¬ 

sionless) Berry flux, with /-^(p) the equilibrium band oc¬ 
cupancy; here we have only kept terms linear in 6n and 
V(/). In addition to the conventional first term, which gov¬ 
erns the behavior of ordinary plasmons, v(r, t) in Eq. (S- 
1^ exhibits a self-induced anomalous velocity component 
Vq that yields chirality as shown in Eig. Note that 
the mass m appearing in Eq. (S-2) is the plasmon mass^ 


which characterizes the collective motion of the Eermi 
sea [23] . 

CBP chirality can be understood intuitively by exam¬ 
ining the anomalous velocity pattern set up by the plas¬ 
mon’s electrostatic potential 0 (a more complete treat¬ 
ment is given below). Due to the cross product in 
Eq. ( |S-2| ), the anomalous velocity flow is directed along 
the equipotential contour lines of (j) (see arrowheads in 
Eigs. and Near an edge, surface charges asso¬ 

ciated with the plasmon wave produce a potential as 
shown in Eig. 0i- The corresponding anomalous velocity 
field directs electrons into the nodal regions to the left 
of each region of negative charge build up (i.e., excess 
electron density), for the orientation shown and iF > 0. 
Thus for a leftward direction of plasmon propagation, the 
anomalous velocity flow assists the collective motion of 
the electronic density wave, leading to faster propagation 
^^dge right-propagating mode the anomalous 

flow works against the collective motion, yielding slower 
propagation 

Crucially, depends directly on the self-generated 

electric field — V(/)(r). Consequently, the magnitude of 
the splitting is governed by the wave vector q and the 
strength of Coulomb interactions. We emphasize, how¬ 
ever, that CBPs are a linear phenomenon, with the mode 
splitting Acd independent of the magnitude of (j). As we 
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FIG. 2 : Chiral Berry plasmons at the boundary of a 
half plane, a) The CBP frequency splitting increases with 
increasing Berry flux, T. Parameter values are as in Fig. 
but with \q\/qQ — 0.3. b) CBP edge-confinement length, 7 ^“ , 
of and see Eq. (|^, for increasing T. For large T^ 

the confinement length of the mode diverges, indicating 

that it joins the bulk, while the mode becomes more 

confined to the edge, (inset) EM profile shown for \T\ = 0.3, 
exhibiting confinement of to the edge. 


will show, the g^-dependent CBP splitting can be large: 

T 1 /S.UJ = ^ A—J^\q\, (3) 

K 


where ^ is a numerical prefactor of order unity that de¬ 
pends on geometry, and we have used the 2D Coulomb 
potential W{q) = 27re/(/^|q|) with background dielectric 
constant k. For edge CBPs, we find A = 8v^7r/9 (see 
below), yielding large splittings hAuj 6 — 60 meV for 
q = 1 — 10/im“^ (here we have used = 1, /i: = D. 
The appearance oi ^ j k on the right hand side of Eq. Q 
signals the crucial role interactions play in Acj. 


Chiral Edge Berry Plasmons - We now analyze col¬ 
lective motion described by Eq. treating the electric 
potential 0 self-consistently. For an infinite bulk, apply¬ 
ing dt to the continuity equation (Ef), using Eq. ( |S-2[ ), 
and substituting in the force equation 0) yields 

/ oo 

dr' W(y — v') dniv'). (4) 

-00 


In arriving at Eq. Q, we have used V -Wa dxdy(j){Y) — 
dydx(t){Y) = 0. Importantly, Berry flux T is absent in 
Eq. and has no effect on bulk plasmon dispersion. 
Indeed, decomposing into Fourier modes 6n ^ r 

and using the Coulomb interaction hF(r, r') = —e/{hi\T — 
r'l) yields the usual 2D plasmon dispersion 


= ao{q/qo), ao = ^^no, qo = (5) 

^ mil 

which remains gapless at g = 0. In contrast, bulk magne- 
toplasmons are gapped due to cyclotron motion [IMT]. 

Close to a boundary, the situation is dramatically al¬ 
tered: Berry flux T leads to the emergence of ehiral 
edge plasmons, one-dimensional chiral analogues of sur¬ 
face plasmons [mEn. This is most easily illustrated for 


an infinite metallic half-plane, where n(r,t) and v(r, t) 
are generically finite for x > 0, but are zero for x < 0. 
Here the plasmon propagates as a plane wave along y, 
with the fields in Eq. 0 taking the form 

0(r, t) = (^q{x)e^^^-^^y, Sn{T, t) = Snq{x)e^^^-^yy. (6) 

The presence of the edge allows charges to accumulate, 
see Fig. [I'd). Indeed, inserting the fields defined above 
into Eq. (^) yields 6n{r)e^^^ = —V • [v{r,t)Q{x)]/iuj. 
Here we have explicitly inserted 0(x) to emphasize the 
vanishing of the velocity density outside the metallic half¬ 
plane [0(x) = 1 for X > 0 and &{x) = 0 for x < 0]. 
Due to the &{x) inside the divergence above, a non¬ 
zero Vx flowing into the boundary induces an oscillat¬ 
ing surface charge component of (5n(r)e*^^ concentrated 
at the edge x = 0 [24] [25]. By replacing Sn{r,t) by 
—V• [v(r, t)Q{x)^/iuj in the integral for 0(r, t) in Eq. (|^, 
we thus find a jump condition for dx(t)q{x) at x = 0: 

9.7|o+-^-<^9Io- = (7) 

where Wq{x) = —{e/n) J^^dk e^^^/\q‘^ + is the 

ID effective Coulomb interaction. Here we decomposed 
Eq. 0 into plane-waves using Eq. 0, and integrated 
across the delta function, dxO{x) = 6{x), which accounts 
for the surface charge layer. The parameter q corresponds 
to the wavevector in y (along the edge), and k describes 
variations in x (perpendicular to the edge). Crucially, 
finite in Eq. ( |S-2| ) makes Va^|o+ depend on both the 
magnitude and sign of q along the edge [45] . 

We now find edge CBP solutions of Eq. 0 using the 
boundary condition 0 and continuity of 0(r). This 
problem is a non-local integro-differential equation owing 
to the kernel Wq{x). We adopt a simplified approach, 
approximating Wq{x) by a similar kernel Wq{x) with 
the same area and second moment [19] |24|: Wq{x) = 
~K f^^dk2qe'^^^l{k‘^ + 2g^). We emphasize that Wq{x) 
is extended in x, and captures the long-range Coulomb 
behavior of Wq{x). Indeed, the Fourier transforms of 
Wq{x) and Wq{x) match for small k/q. This method has 
been used successfully to mimic Coulomb interactions in 
isolated systems [T6j [26] . 

Crucially, the simple form of Wq{x) above allows the 
integro-differential equation in Eq. Q to be expressed as 
a purely differential one with 0g(x) obeying 

(.dl - 2q^)(j)>ix) = ^Aslsng{x), {d^ - 2q^)(f)<{x) = 0, 

Kj 

(8) 

where (j)^ and (j)^ are defined inside {x > 0) and outside 
{x < 0) the sample, respectively. Eq. ^ yields simple <p 
profiles 


(j><{x) = 


(f>> {x) =(f>ie 


(9) 
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where 70 = \/ 2 |g|, and 71 = V2\q\ - 

co’^]/[2(co’q^^^)^ — The latter was obtained from 

Eqs. Q and ® by eliminating Suq. 

Using Eq. the boundary conditions of continuous 
(j) and the jump condition 0 can be expressed compactly 
via the relation 54> = 0 , with 4> = (^o, 0 i)^- 


Eig. 1^ inset), suppressing its coupling to bulk plas- 
mons. Therefore, above the threshold ujth (gray region in 
Eig. [^), the single, well-defined, mode propagates 

unidirectionally along the edge. When hybridized with 
light, it will allow for strong non-reciprocal propagation 
of CBP-polaritons without magnetic field (see below). 


5 = 


-1 


\/2\q\ 7 i - -D 


D = 


271 ( 7 "“^)" q^J^sgn{q) 




UJ 


( 10 ) 


where D was obtained using Eqs. (S- 2 ) and Q, and 
F = j kTi. Left- and right- moving plane wave 

modes along the edge, can be identified through 

the zero modes of S. We first note that for = 0 , solv¬ 
ing det(5) = 0 yields degenerate non-chiral edge modes 
with = (2/3)see Refs. pTl [T^ [24) . 

Eor non-zero T, the zero mode solutions of 


Eq. ( [Tq| ) become split, yielding chiral edge plasmons 
(CBPs) as shown in Eigs.lfjand^ The modes and 
^^dge wuvcs in opposite directions along the 

edge, with faster and slower speeds, respectively. Impor¬ 
tantly, the frequencies depend both on q along the 

edge (Eig. 0 ), as well as T (Eig. [^). The splitting be¬ 
tween mo des g rows with q and T^ since = eV(l)X F/Ti 
[see Eq. ( |S- 2 | ]. Indeed, for small qF (so that qF 
^buik), obtain an approximate dispersion for edge 
CBPs as ( 2 / 3 )i/ 2 ,^buik ^ y/2\q\F + 0{q^F^). 

As a result, we obtain the g^-dependent Auj in Eq. 0 . 
This behavior sharply contrasts with that of magneto- 
plasmons, which have a g^-independent splitting given by 
the cyclotron frequency [n m na m]. As a result, 
far larger splittings, arising from interactions, can be 
achieved for CBPs. 

Interestingly, for large enough q and/or F^ the 
mode (blue line in Eigs. [^,b and [^) merges with the 
bulk plasmon mode (dashed line). As this mode 

merges with the bulk, its potential profile ceases to be 
localized along the edge. This is shown by a diverg¬ 
ing confinement length of the electric potential, 7 C ^5 see 
Eig.j^. In contrast, stays far from the bulk disper¬ 

sion, yielding a potential (and electric field) tightly con¬ 
fined to the edge. The threshold cjth above which 


merges with the bulk can be obtained from Eq. (10). Set¬ 
ting UJ = in Eq. (10) yields the threshold frequency 


hujth = 


firiQ 

V2m\F\ 


^12 


= 18.3 


no [cm ^]/10 
(m[me]/0.03) x \F\ 


meV, ( 11 ) 


with rrie the free electron mass. Eor scale we consider a 
plasmon mass 0.03me, as measured in graphene [23] . 

Conservation of uj and q along the edge protect the 
^^dge from coupling to bulk 2 D plasmons. Scat¬ 

tering processes that relax q contribute to propagation 
losses. However, the tight edge confinement of the uj^^^ 
mode suppresses its electric field in the bulk regions (see 


Experimental Signatures of CBPs - Strong plas¬ 
mon mediated light-matter interactions [IIHlS] make op¬ 
tics an ideal means of probing/controlling CBPs. Pho¬ 
ton coupling to plasmons with gapless dispersion (e.g., 
2D plasmons, and CBPs here) can be achieved through 
strategies such as gratings, and prism geometries HB. 
Observing unidirectional (non-reciprocal) propagation in 
such setups can reveal the existence of CBPs. Eor demon¬ 
stration, we detail an alternative experimental probe: 
CBP-photon coupling in finite geometries, such as metal¬ 
lic disks, where dipolar plasmonic modes can dominate 
optical absorption muni. 

In metallic disks with finite F, CBPs manifest as 
clockwise/anti-clockwise moving plasmonic dipole modes 
(Fig. § 1 ). These modes can be described via a sim¬ 
ple oscillator model for the motion of the dipolar CBP 
center of mass (COM) coordinate [22], {x}, where {•} 
denotes the COM average. Here {v^} ~ F {eV0} x z 
(green arrow) gives rise to an intrinsic angular frequency 
uja of plasmons in a disk (orange arrow), which adds 
to (subtracts from) the plasmon frequency ujq to pro¬ 
duce non-degenerate anti-clockwise (clockwise) rotating 
modes (Eig. [^); see detailed derivation in SOI. Here we 
have used F pointing to positive z. A bosonic analogue 
for ultra cold atomic gases is discussed in Ref. m- 

With an a.c. probing electric field ^ the COM 

equations of motion are: 

d‘^{x} + LO^ {x} + C0adt{y} = -eEx 

dt{y}-‘^adtix] 110% 

( 12 ) 

Here ujQ{d) is the bare plasmon frequency in a disk of 
diameter d, in the absence of Berry curvature. 

Writing the current density as j = eng 9t{x}, we invert 
the COM equations of motion to obtain the optical ab¬ 
sorption (real part of the longitudinal conductivity [22]). 
As shown in Eig. [^, we find a split peak structure with 
the dipolar CBP peaks given by the poles of Eq. ( [T^ : 

^disk ^ = FTiuj^ « meV, 

V 4 2 d[/imj 

(13) 

where Auj = and Tiuj^ = uj^m/riQ. 

On the right side, we have estimated ujq ~ 27re^no|q|/m, 
with |q| 1/d (approximating the lowest lying plasmonic 

excitation in a disk) [46]. Here we have used k = 1. 
Importantly, Auj depends on the disk size, d, a unique 
property of CBPs. 

The tunable optical absorption split peak structure 
(via F and d) in the absence of an applied magnetic 
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field gives a clear experimental signature of CBPs. In 
plotting Fig. Ht, we have included the damping rate phe¬ 
nomenologically via ^ + FQt, yielding a Lorentzian 
lineshape with its half-width determined by F. Split 
peaks are clearly visible when Auj ^ F, yielding peaks 
at To give a sense of scale, we note a typical value 

hr ^ few meV, see e.g. Ref. [28] where hV ^ 4 meV was 
measured in graphene disks. Using F/cj* = 0.25 and tak¬ 
ing = 1 , clearly resolved peaks can be resolved 

for disk sizes d ^ 1 /im (Fig. ^). 

CBPs in Anomalous Hall Materials - We now dis¬ 
cuss materials where CBPs can be realized. We predict 
that metallic systems with non-vanishing will generi- 
cally host CBPs. Finite T requires broken time reversal 
symmetry, and may arise in magnetically ordered systems 
or out-of-equilibrium non-magnetic systems (see below). 
The former includes magnetically doped semiconducting 
quantum wells (see e.g.. Ref. m, where ^ 1/2 was 
predicted) and topological insulators jSHSj. 

As a concrete example, we examine the magnetically 
doped topological insulator chromium doped thin-film 
(Bi,Sb) 2 Te 3 , which was recently experimentally real¬ 
ized [5]. When moderately doped with electrons or holes, 
we estimate that T ^ 1 can be achieved based on the 
measured anomalous Hall conductivity ^ [5|. 

This yields a large splitting hAuj « 10 — 100 meV for 
short wavelength plasmons, see Eqs. © and When 
probed in the disk geometry, we predict this system will 
yield two split absorption peaks in the absence of mag¬ 
netic field. 

“On-demand” CBPs in Non-Magnetic Materi¬ 
als - Intriguingly, finite can also be achieved in non¬ 
magnetic materials, without an applied magnetic field. 
This includes, for example, gapped Dirac materials where 
inversion symmetry is broken. In these, a valley polariza¬ 
tion (Fig.[^ inset) can be induced by circularly polarized 
light [30], yielding an anomalous Hall effect that has re¬ 
cently been observed [7|; see materials discussion below. 

Can CBPs exist in photoexcited systems? To ana¬ 
lyze this, we focus on nominally time reversal invari¬ 
ant gapped Dirac materials. We model the valley de- 
pendent as [H]; J^k,k' = sgnA, 

where uk^k' are the valley carrier densities, Tz = T for 
the A, K' valleys, = 2 is the spin degeneracy, and 
h = A?/4'Kv‘^h^ gives a characteristic density scale. The 
bandgap is 2 |A|. When uk = nx'^ as in equilibrium, the 
total flux T = Tk + ^'k vanishes. 

Interestingly, when the system is pushed out of equi¬ 
librium, e.g., by circularly polarized light, the popula¬ 
tions in the valleys may become imbalanced, uk tik' 
122130] (Fig.|^ inset). As a result, the net Berry flux for 
the entire electronic system (summed over both valleys) 
is non-zero, yielding 7 ^ 0. We analyze the collec¬ 
tive modes of the photoexcited electron-hole system in 
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FIG. 3: Chiral Berry plasmons in a disk and in valley 
polarized gapped Dirac systems, a) Illustration of CBPs 
in a disk, showing anti-clockwise/clockwise rotating dipole 
modes. The mode splitting arises from the intrinsic angu¬ 
lar frequency, cua, induced by the combination of T and the 
self-induced electric field, b) Optical absorption split peaks 
for d ipo lar CBPs in a disk (panel a) obtained by inverting 
Eq. ( [T^ for T — 0,0.5,1.0 (dashed black, purple, orange). 
Smaller disk sizes and/or larger T produce larger splittings, 
see Eq. (13). Parameter values: co^o/cu* = 1 and Vju:^ — 0.25. 
c) CBP dipole modes for valley polarized gapped Dirac sys¬ 
tems, with nK — n.pe, and n^' — 0. Characteristic density is 
h — A"^/Aivvph^. (inset) Valley polarization may be induced 
by above-gap circularly polarized light. 


the disk geometry following Eq. (12) above, accounting 
for the mutually attracting electron and hole populations 
and using T = This yields two chiral CBP modes, 
= + 2 co’q ± cja, see SOI. Estimating cjo as be- 

curves in Eig. 3b with 


low Eq. (13), we obtain the 

given by Eq. (13). The frequencies of both modes van¬ 
ish at zero photoexcited carrier density, ripe- The mode 

and can reach sizable values 


^ ,el-h 


,el-h 


splitting increases with Upe, 
for Upe > n, reflecting J^pe 1 
Large splittings require n 


> 
pe rsj 


h. The large gaps 


A ^ leV of many transition metal dichalcogenides (T] 
yield large characteristic densities requirements. 


^ 10 ^^ cm In contrast, other gapped Dirac ma¬ 


terials such as G/hBN [34ll36] and dual-gated bilayer 
graphene m possess A 10 — 200 meV, yielding sig¬ 
nificantly smaller and more favorable h. These materials 
possess valley-selective optical selection rules [22l|30], and 
present an ideal venue to achieve maximal T = 1 (and 
large Auj), even with relatively weak pump fluence [47] . 

A further promising strategy to achieve large CBP 
splittings is to stack m layers of gapped Dirac materi¬ 
als on top of each other, with no tunnel coupling be¬ 
tween the layers. Stacking achieves (i) larger photo ex- 
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cited carrier densities due to the increased absorption, 
and (ii) a larger maximal Berry flux (when Upe ^ ft) and, 
hence, larger CBP splittings; the maximal Berry flux is 
‘^max = ^‘^singie* such a structure, the long-range 
Coulomb interaction allows the photoexcited densities in 
different layers to oscillate collectively. 

Conclusion - CBPs are robust collective excitations 
of metallic systems, arising from two simple ingredients: 
Berry flux and interactions. Our analysis indicates that 
CBPs survive for both weak and strong interactions. As 
a result, we conclude that CBPs are generic in metal¬ 
lic anomalous Hall phases, including out-of-equilibrium 
states where finite emerges from driving (e.g., in opti¬ 
cally pumped valley polarized gapped Dirac materials). 
Indeed, this feature allows CBPs to be used for all-optical 
diagnostics of anomalous Hall and topological phases, as 
well as pumped or periodically driven systems, e.g., as in 
Ref. [381 EHl • Optical probes of the latter are particularly 
appealing since transport measurements require contacts 
that often complicate and destroy the novel electron dis¬ 
tributions and coherences of driven systems. 

Perhaps the most appealing prospect is coupling CBPs 
with light for sub-wavelength, and strong non-reciprocal 
propagation. Exhibiting a single chiral mode at large 
q (large splitting, and large frequency), precisely where 
plasmons give large compression of optical mode volume, 
hybrid CBP-polaritons are strongly non-reciprocal m- 
As we propose, CBP mediated unidirectional waveguides 
can be realized in readily available non-magnetic materi¬ 
als (e.g., the van der Waals material class). CBP based 
non-reciprocity, if realized experimentally, stands to play 
a vital role in the miniaturization of optical components 
that are magnetic field-free. 

Note added: During the review of our manuscript, a 
related work on CBPs in photoexcited gapped Dirac sys¬ 
tems appeared, see Ref. m- 
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SUPPLEMENTARY ONLINE INFORMATION 
FOR “CHIRAL PLASMONS WITHOUT 
MAGNETIC FIELD” 


In this supplement we include discussions of an alter¬ 
native formulation of chiral Berry plasmons based on cur¬ 
rent densities and the conductivity tensor, chiral Berry 
plasmon dipole modes in a disk, optical selection rules for 
gapped Dirac materials with inversion symmetry break¬ 
ing, and photo-induced chiral Berry plasmons. 


Current density and the conductivity tensor 


For a complementary view to the approach employed 
in the main text, here we present an alternative formu¬ 
lation of the plasmon equations of motion. Rather than 
focusing on the semiclassical equations of motion for the 
particle and velocity density fields, here we work with 
the charge density p(r, t) = — en(r, t), the current density 
j(r, t), and the conductivity tensor a. As we will be in¬ 
terested in oscillatory collective (plasmon) modes, it will 
be useful to work in terms of Fourier modes, p(r, t) = 

and = </><^(r)e*“‘. 

The charge and current densities obey a continuity equa¬ 
tion, which in terms of the Fourier modes reads: 


icjp^(r)-|-V-jc^(r) = 0, j^(r) = t7(cj)V(/)c^(r). (S-1) 


Importantly, even at 5 = 0, the conductivity a{uj) pos¬ 
sesses off-diagonal contributions, captured by the anoma¬ 
lous Hall conductivity 

In terms of frequency dependence, the longitudinal 
conductivity takes a standard Drude form, while the 
Hall conductivity, arising from the anomalous velocity 
X ri(p) [8], is frequency-independent: 




rioe^ jm 


lUJ 


A 

CT™,, = const. 


(S-2) 


Here uj = uj — where Ttr is a transport scattering 

time that arises, e.g., from impurity scattering. We con¬ 
sider collective modes in the regime ooTtr ^ 1, where the 
1/rtr contribution to axx can be neglected. 

The anomalous Hall conductivity is easily connected 
with quantities in the main text via jh. Here 

T = f <^^P^i(p)/i^(p)/(27rfi)^ is the dimensionless 
Berry flux, with f^(p) the equilibrium band occupancy, 
and ri(p) the Berry curvature.Throughout this work we 
only consider terms up to linear order in V(/). 

In this formulation, chiral Berry plasmons (CBPs) at 
system edges or in disks arise in exactly the same way 
as shown in the main text, with v(r, t) and T tra ded for 
j(r,t) and That is, Eqs. (S-1) and (S-2) of this 

supplement, with 1/rtr ^ 0, are equivalent to Eqs. (1) 
and (2) of the main text; following the same steps of 
derivation yields Eqs. (4-11) of the main text with the 
replacement T = 


We note parenthetically that Eq. (S-2) is different from 
what would be obtained in the analogous case for magne- 
toplasmons in a system subjected to an applied magnetic 
field. While both cases feature a nonzero Hall conduc¬ 
tivity (Txy^ in the case of magnetoplasmons the applied 
magnetic field induces a Lorentz force that affects elec¬ 
tronic motion by modifying the force balance. This leads 
to different forms of cFxx and (Jxy from those in Eq. (S-2); 
in particular, (Jxy in the presence of a finite magnetic fieid 
exhibits an uj dependence [24], in contrast to Eq. (S-2). 
This yieids quaiitativeiy different behavior of coiiective 
modes. Eor exampie, buik magnetopiasmons are gapped, 
whereas buik piasmons in the presence of a Berry fiux 
giving the same dc-Haii conductivity are gapiess (see dis¬ 
cussion in main text). 


CBP dipole mode in a disk 


Here we present a more complete analysis of the CBP 
dipole modes in a disk, which dominate the optical ab¬ 
sorption in metallic disks. These dipole modes can 
be conveniently described through the center of mass 
(COM) motion, {x(t)}, wherein all internal forces cancel 
(viz. Newton’s third law). Here {•} denotes the COM 
average, with {x(t)} = f d^xn(x,t)x, and {p(t)} = 
f d^xp(x,t). The COM equations of motion can be ob¬ 
tained from Eq. (1) of the main text via integration by 
parts, along with the condition that the velocity vanishes 
when |x| ^ oo. This yields 


at{x} = 1^1 + ^{eV<p X f2(p)}, at{p} = {eV(/)}, 

^ (S-3) 

where E is a self-generated electric field associated with 
the plasmon motion. 

To capture the dipolar mode, we use the harmonic po¬ 
tential —e(j){x) = electrons (see e.g.. Ref. 

[28]) where cjq is the bare plasmon frequency in a disk 
(with diameter d) in the absence of Berry curvature. 

In analyzing the COM equations of motion, we first 
note that eV0(x) depends only on x, and ^4(p) only on p. 
Keeping terms to linear order in 5n^ we find the anoma¬ 
lous velocity contribution {v^} for the COM as 


/ ^2 




X z C{x(t)} X z, 

(S-4) 


where the Berry flux is given by J^(x, t) = 
^ / d^p Ui(p)/i(x, p, t)/( 27 rfi)^, the constant C can be 
obtained as detailed below, and i is a band index. Here 
we will concentrate on the motion of bulk electrons in 
a given band, n. We therefore take /i<n = 1 every¬ 
where inside the disk for bands below n (but vanishing 
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outside the disk). For the purpose of estimating param¬ 
eters, we make an assumption of local equilibrium and 
set /n(x, p,t) = + 1]“^, where /i(x, t) is a 

space and time varying chemical potential. Adopting a 
simple model of a rigidly moving disk of charge with con¬ 
stant density no, and density equal to zero outside, gives 
T = jriQ and Here T is obtained 

with a fixed uniform chemical poten tial j i. 

Using {eV0} = —Eq. (S-4), and substitut¬ 
ing into Eq. (|S-3|), we obtain the equations of motion for 


the dipole mode: (5| + Aij){xj} = 0, where 


A=( ‘"o' ‘"“f* 

-UaOt CJg 




no?i 


(S-5) 


Here we have used + kept only 

terms up to linear order in Sn. 

Writing {x} = xoe*^^, we obtain a secular equation 
Mxq = 0 where M = —uj‘^ + A, with dt replaced by 
iuj within A. Plasmons are given by the zero modes, 
det(M) = 0, yielding the split dispersion relation 


where uj±, ujq > 0. (In the following we will use only posi¬ 
tive frequencies; a similar analysis yields the same results 
for negative uj±,ujo as well). These plasmon modes are 
chiral as they correspond to a rotating COM displace¬ 
ment and momentum 

Mt)U = ^ ( ii) 

(S-7) 

where {p(t)}± and {x(t)}± are offset by a phase of 7rl2. 
As a result, a rotating {x(t)}± gives rise to a circulat¬ 
ing momentum/current density (see Fig. 3a of the main 
text). The clockwise and anticlockwise (a;+ and uj-) mo¬ 
tions of chiral plasmons sketched in Fig. 3a are oriented 
for a positive Berry flux, pointing in the z direction. 
This direction sets the orientation of uja- For opposite 
sign of iF the orientations are switched. The distinct 
frequencies uj± arise from the anomalous velocity {v^} 
adding/subtracting propagation speed from the modes 
at iF = 0. 

The CBP dipole mode can manifest in distinct split 
peaks for optical absorption. This splitting can be ana¬ 
lyzed by writing the current density as j = eno9t{x} by 
inverting M above, and relating current to electric field 
via j = gE, where E is probing field, and g is the con¬ 
ductivity tensor. Optical absorption is characterized by 
the real part of the longitudinal conductivity, gxx, as 

= - ^ (a;2 ± ’ (§-8) 


2 

where V = is the Drude weight, and F is the 

transport relaxation rate, included phenomenologically 


via ^ FTdf. This yields the split peak optical 
absorption for the disk geometry shown in Fig. 3b of the 
main text. 


Optical Selection rules for gapped Dirac materials 

Here we detail how circularly polarized light can yield 
valley selectivity in gapped Dirac materials (GDMs) 
where inversion symmetry has been broken. This selec¬ 
tivity arises due to the pseudo-spinor nature of the wave- 
functions on the A/B sublattices. The selection rules 
were derived in Ref. m by analyzing the orbital mag¬ 
netic moments of electronic wavepackets in the K and 
K' valleys. For the reader’s convenience, here we present 
an alternative calculation based on Fermi’s golden rule. 

The low energy Hamiltonian for gapped Dirac materi¬ 
als can be described as = Hk + Hk' , where 


Hk = VCTj^ • Pk + ^CFz, Hk' = VCT- • pk/ + Acr;^, 

(S-9) 

where Hk^ Hk' describe electrons close to the K and K' 
points, V is the Fermi velocity, and cr± = aa^xTcr^y, with 
Ike Pauli matrices. Here y^k,k' describe momenta 
taken relative to the K and the K' points, respectively, 
and 2A is the gap size. In the following we shall drop 
the explicit K and K' labels on Pk,K' for brevity. 

A variety of systems obey Eq. ( |S-9[ ) , including van der 
Waals heterostructures where A/B sublattice symmetry 
has been broken as in G/hBN heterostructures. In such 
systems the gap size 2A corresponds to the asymmetry 
of the potential on the A/B sublattices. Below we will 
focus on this case for concreteness. However, the un¬ 
derlying physics is general and applies to a broad range 
of GDMs, such as dual-gated bilayer graphene and tran¬ 
sition metal dichalcogenides, where inversion symmetry 
has been broken. 

The eigenfunctions of Eq. (S-9) can be expressed as 
pseudo-spinors 


sinf \ 
—cosf J ’ 
(S-10) 

where |±) denote states in the conduction (valence) band, 
r = 1 for valley K and r = — 1 for valley A', tan 0 = 
u|p|/A, and tan (/) = py/p^. Here the energy eigenvalues 
are . 

We proceed by noting that when light with frequency 
huj > 2A is incident on the GDM, electrons in the valence 
band can be excited into the conduction band. The light- 
matter coupling is captured by writing p p — eA/c, 
where the vector potential A is related to the incident 
light electric field E via A = ^E. Here E = Eqc*^^ for 
light with frequency uj. The rate of electron-hole pair 
creation (absorption of photons) can be calculated via 




cos|e-*^'^ 

• 

sin 2 


\-)k{K') = 
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Fermi’s golden rule 

WKiK') = Y E - nu/2), (S-11) 

P 

with the matrix elements 


However, a non-vanishing net Berry flux for the 
photoexcited system is achieved when uk 7 ^ tik' (Fig. 3c 
inset of the main text). Setting uk' = 0, giving Tk' = 
+1, the valley polarized electron (nei) and hole (nh) pop¬ 
ulations concentrated at the K valley conduction and va¬ 
lence band extrema yield 


= —{+\E,a, + Eyay\-)K 

UJ 

= —{+\E,cj, - Eyay\-)K'. (S-12) 

UJ 


We note that the different signs in front of Ey in and 
arise from the different way A couples to pseudospin 
in Hk and Hk' (see Eq. S-9). 

With the help of the identities 


{+Wx\-)k(K') = single - cos^f(S-13) 


and writing Eq = |Eo| (x + i? 7 y )/\/2 for left-handed (LH, 
7 ^ = 1) and right-handed (RH, 7 ^ = — 1) circularly polar¬ 
ized light, we obtain 


Wl = Wo 


■2A 

2 

'2A 


, W^^, = Wo 



(S-14) 


Here Wo = Ase^|Eop/(16?i^co’), where Ns is the spin de¬ 
generacy. 

Equation ( |S-14[ ) shows that electron-hole transitions in 
the valleys KjK' can be selectively excited using LH/RH 
circularly polarized light (see Eig. 3c inset of the main 
text). Indeed for Tiuj = 2A, perfect selection of K or K' 
electron-hole transitions can be achieved, in agreement 
with Ref. [30]. We emphasize that this selectivity comes 
from the orbital physics of light-matter coupling; it does 
not require or involve spin-orbit coupling and can even 
arise in materials with negligible spin orbit-coupling as 
modeled by Eq. (S-9). 


Optically pumped “On-Demand” CBP dipole mode 
in gapped Dirac materials 

Here we consider how CBPs might arise in non¬ 
magnetic gapped Dirac materials such as M 0 S 2 , 01 * 
gapped G/h-BN. To analyze this, we focus on nomi¬ 
nally time r ever sally invariant gapped Dirac materials. 
We model the valley dependent T as [31]: Tk,k' — 

~ 1/2 

sgnA, where uk.k' are the valley car¬ 
rier densities, = T for the A, K' valleys, = 2 is the 
spin degeneracy, and n = Hjj ?gives a characteris¬ 
tic density scale. The bandgap is 2|A|. When uk — 
the total flux E = Ek + E'j^ vanishes. 




pe 


1-1 


11/2 


il/2 


+ 


n = 1.8 X 10 


2 ly/n + 77h ^/n + 77el- 

13 (A[eV])^ 

(7;i?[cms“^]/10^)^ 


cm 


(S-15) 


Here we have used that Berry curvature (including its 
sign) is the same for electrons and holes, following from 
particle-hole symmetry of the Dirac Hamiltonian, and 
used the convention 77ei,h > 0. In deriving Eq. (S-15), we 
used the spin degeneracy Ns = 2 for the K' valley, and 
noted that only a single spin species in the K valley is 
excited by the circularly polarized light. 


Eor demonstration we examine CBPs in the disk ge¬ 
ometry as above. We analyze the coupled motion of pho¬ 
toexcited electrons and holes in a single valley via COM 
coordinates ({xei}) and ({xh}), respectively, giving: 




(S-16) 


where Aei and Ah are defined as in Eq. (7) of the main 
text, with the density no and plasmon mass m replaced 
by the appropriate values for electrons or holes. In writ¬ 
ing these equations, we have used that the restoring force 
arising from the mutual attraction of photoexcited elec¬ 
tron and hole subsystems is {eV0}ei(h) ~ T<^({xei} — 
{xh}), where {•}ei,h denote the COM averages for elec¬ 
tron and hole distributions and the upper (lower) sign is 
for electrons (holes). Here a characterizes the strength 
of the electron-hole interaction. Eor brevity, in the fol¬ 
lowing analysis we set Aei = Ah, giving a = mujQ, where 
ujQ is the plasmon frequency associated with a unipolar 
system with carrier density n = riei = ^h- 


Zero modes of Eq. (S-16) with 
with {xei} = -{xh}, giving: 


yield CBPs 


^ ^a 2cJq ± UJa, (S-17) 


where uja is given by Eq. (7) of the main text with E re¬ 
placed by Importantly, chiral (electron-hole) 

plasmons feature co-rotating {xei} and {xh}. This dipole 
like rotation follows from the equal sign of E for electrons 
and holes. Modes with {xei} = {xh} are also eigenmodes 
of Eq. (S-16). However, they have frequency, u; = 0. 
























